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We show that spin generalization of elliptic Calogero-Moser system, elliptic extension 
of Gaudin model and their cousins can be treated as a degenerations of Hitchin systems. 
Applications to the constructions of integrals of motion, angle-action variables and quan- 
tum systems are discussed. 

The constructions are motivated by the Conformal Field Theory, and their quantum 
counterpart can be treated as a degeneration of the critical level Knizhnik-Zamolodchikov- 
Bernard equations. 
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1. Introduction 

Integrable many-body systems attract attention for the following reasons: they are 
important in condensed matter physics and they appear quite often in two dimensional 
gauge theories as well as in conformal field theory. Recently they have been recognized in 
four dimensional gauge theories. 

Among these systems the following ones will be of special interest for us: 

1 . Spin generalization of Elliptic Calogero - Moser model - it describes the system of 
particles in one (complex) dimension, interacting through the pair-wise potential. The 
explicit form of the Hamiltonian is: 

N 2 

i=\ i^j 

where Zi are the positions of the particles, pi - corresponding momenta and Si are the 
"spins" - / x / matrices, acting in some auxilliary space. The conditions on Si will be 
specified later. The only point to be mentioned is that the Poisson brackets between 
p, z, S are the following: 

{Pi,zj} = 8ij 

{{Si)ab, (Sj) c d} = d~ij(d~ad(Si)bc — 0~bc(Si) a d) 

2 . Gaudin model and its elliptic counterpart. We describe first the rational case. 
Consider a collection of L points on P 1 in generic position: w±, . . . , wl, assign to each wi 
a spin Si (an N x N matrix) and define the Hamiltonians [G] : 

Tt(SiSj) 



The main goal of this note is to include these two (seemingly) unrelated models in the 
universal family of integrable models, naturally related to the moduli spaces of holomorphic 
bundles over the curves. It will turn out, that the appropriate objects to study are Hitchin 
systems. As a by-product we shall invent elliptic Gaudin model, which includes both cases 
as a special limits. We shall also obtain a prescription for construction of integarls of 
motion and action-angle variables. The paper is organized as follows. In the section 2. we 
remind the construction of Hitchin systems. The section 3. is devoted to the explanation 
of the mapping between the Hitchin systems and the models, just described. The section 
4. deals with action-angle variables and integrals of motion. We conclude with the remarks 
on the quantization of our constructions. 
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3. Construction of the Systems 

We must confess that all the models we are discussing are motivated by the studies 
of Knizhnik-Zamolodchikov-Bernard equations [KZ],[Be], [Lo] , [F V] , [I] . Our paper is a de- 
velopment of [NG1], [NG2] . One of the outcomes of our work might be an insight in the 
W- generalizations of them. 

3.1. Hitchin systems 

Hitchin has introduced in [H] a family of integrable systems. The phase space of 
these systems can be identified with the cotangent bundle T*M to the moduli space M of 
stable holomorphic vector bundles of rank N (for GL N (€) case) over the compact smooth 
Riemann surface E of genus g > 1. His construction can be briefly described as follows. 
Fix the topological class of the bundles (i.e. let us consider the bundles £ with c\{£) = k, 
with k - fixed). Consider the space A s of stable complex structures in a given smooth 
vector bundle V, whose fiber is isomorphic to € N . The notion of stable bundle comes from 
geometric invariants theory and implies in this context, that for any proper subbundle U: 

degjU) deg{V) 
rk(U) rk{V) 

The quotient of A s /G of the space of all stable complex structures by the gauge group is 
the moduli space M. Its dimension is given by the Riemann-Roch theorem 

dim{N) = N 2 (g-1) + 1 

Now consider a cotangent bundle to A s . It is the space of pairs: 

(f), d" A 

where is a Matjv(C) valued (1,0) - differential on E, d" A is an operator, defining the 
complex structure on V: 

d" A :Q°(E,V)^Q°' 1 (E,V) 
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The field <p is called a Higgs field and the pair d A , 4> defines what is called a Higgs bundle. 
In the framework of conformal field theory the Higgs field is usually reffered to as the 
holomorphic current, while holomorphic bundle defines a background gauge field. 

The cotangent bundle T*A S can be endowed with a holomorphic symplectic form: 

Tr5(/> A 5d A 

where 5d A can be identified with a (0, 1) - form with values in N by N matrices. Gauge 
group Q act on T*A S by the transformations: 

-> 9~ X <\>9 
d" A -> g-'d'^g 

and preserves the form u. Therefore, a moment map is defined: 

H = [d% <f>] 

Taking the zero level of the moment map and factorizing it along the orbits of Q we get 
the symplectic quotient, which can be identified with T*M. Now the Hitchin Hamiltonians 
are constructed with the help of holomorphic (1 — j, 1)- differentials Uj^. where ij labels a 
basis in the linear space 

H 1 ^^®^) =C (2j '- 1)(9_1) 

for j > 1 and <C 9 for j = 1. Take a gauge invariant (j, 0) -differential Tr^ and integrate it 
over E with the weight Vjj. : 

Hj,i 3 = [ ^,Tro' 
is 

Obviously, on T*A S these functions Poisson-commute. Since they are gauge invariant, 
they will Poisson-commute after reduction. Also it is obvious, that they are functionally 
independent and their total number is equal to 

N 

g + ^(2j - 1)(<7 - 1) = iV 2 G7 - 1) + 1 = dim(Af) 
Therefore, we have an integrable system. 
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3.2. Holomorphic bundles over degenerate curves 

Now let us consider a degeneration of the curve. Recall, that the normalization of the 
stable curve E is a collection of a smooth curves E a with possible marked points, such 
that any component of genus zero has at least three marked points and every component 
of genus one has at least one such a point. For each component E a we have a subset 
X a = {x a , . . . , x^ a } of points. Let us denote the pair (E a , X a ) as C a . The disjoint union 
of CaS is mapped onto E by the normalization map n. Let us denote by X a p the set of 
double points 7r(X a ) C\ir(Xp) for a ^ f3 and as X aa the set of double points in 7r(X Q ) (these 
appear due to pinching the handles). The union of all X a p we shall denote by X C E. 
We define x l £p £ X a p as ir(x l a ) fl 7r(x^). Notice, that it may be empty. Stable bundle S 
over E is a collection of holomorphic bundles £ a over E a of rank N (there might be some 
generalizations with different ranks of the bundle over different components - these are 
unnatural as a degeneration of the bundle over smooth curve) with the identifications g^p 
of the fibers 

with the obvious condition: g l ^pgp l a = 1. 

Gauge group acts on the complex structure of the bundle £ a for each a as in the 
smooth curve case. The new ingredient is the action on g 1 ^. Fix a gauge transformations 
g a for each component of E. Then g^p are acted on by g a as follows: 

fa/3 ^(^) _1 fla>a(^) 

Now we have to introduce a notion of stable bundle. The condition of stabilty is: 

For each collection of proper subbundles T a C £ a , such that 

and 

rk{F a ) = N' < N 
for each a the following inequality holds: 

N' 

deg{T a ) < —deg(£ a ) 

for any a. 
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Let A will denote space of collections of d!' A a operators in each £ a together with gr^ 
for each a and f3. Let A s will denote the subspace of A, consisting of stable objects. The 
cotangent bundle T*A S can be identified with the space of collections of pairs 

(£ a , </> a ), <j> a e H 1 ' ^) <g> End{£ a ) 

and 

(9a^Pa,^PaP G T > Hom{£ a |^ , SpU ) 

We normalize = ~Ad*(g l £p)p 3 p a . The Higgs fields a are allowed to have sin- 

gularities at the marked points. As we will see, they could have poles there. Now we 
shall proceed as in the previous section. Consider the gauge group action on T*A S . Since 
the gauge group Q is essentially the product of gauge groups Q a , the moment map is a 
collection of the moment maps for each component E a : 

^ = [^^ + 1?^) 

where the sum over i runs from 1 up to L a while (3 and j are determined from the condition, 
that rr(x J f3 ) = 7t(x l a ). Let us now repeat the procedure of reduction. At the first step we 
should restrict ourself onto the zero level of the moment map. It means, that 4> a becomes 
a meromorphic section of the bundle End(£ a ) £§> Q 1,0 (T, a ) with the first order poles at the 
double points. The residue of <p a at the point x l a equals to for appropriate This 
condition is compartible with the definition of the canonical bundle over the stable curve. 
On the next step we take a quotient with respect to the gauge group action and get the 
reduced space T*Af. The space Af is the quotient of A s by Q. The symplectic form on 
T*Af can be written as: 

a (a,i),(P,j) 

Let us calculate the dimension of T*J\f. We shall calculate the (complex) dimension of M 
by means of the following trick. The moduli space M can be projected onto the direct 
product of moduli spaces Af a of the stable bundles over E Q 's. Actually, the map is to the 
product of the moduli of holomorphic bundles, but the open dense subset, consisting of the 
stable bundles is covered. The projection simply takes the collection of £ Q 's to the product 
of equivalence classes in A/" Q 's. The fiber of this map can be identified with the quotient 
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G/H, where G is the set of collections of g^p, while H is the group of automorphismes of 
x a £ a . This group is a product over all components E Q of genus g(T, a ) < 2 zero of the 
groups H a . For genus zero component H a is GLat(C), while genus one component provides 
a maximal torus - (<C*) N . Therefore, at generic point, we conclude, that the dimension of 
N is 



dim{N) = dim(N a ) + dim(G/H) = N 2 E(Z) + ^ N 2 (g(Y, a ) - 1) = N 2 (g - 1) + 1 



E(L) is the total number of double points. 

3.3. Hamiltonian systems on T*Af 

Now we shall define the Hamiltonians. For each a we take v a ,l,k - the /c'th holomorphic 
(1 — /, 1) differential on E a — X a and construct a holomorphic function on T*A S : 



Obviously, all H aj i : k descend to T*M and Poisson-commute there. 

One can notice, that the integrable systems we get can be restricted onto the invariant 
submanifolds. Namely, the conjugacy classes of are invariant under the flows, which 
we have constructed. Indeed, the casimirs Tr(p^)' are the coefficients at the most singular 
terms of Tr<^, therefore, they are among the Hamiltonians we have constructed. 

4. Gaudin model, Spin Elliptic Calogero-Moser System and so on ... 

4-1. Genus zero models 

Consider a component of genus zero. Let us describe explicitely the part of T*Af 
related to this component as well as the Hamiltonians. We shall omit the label a 
in the subsequent formulas to save the print. Thus, we have C = (P , x\, . . . , xl)- 

Gaudin model 

Assume first, that for i ^ j: Tr(xi) ^ n(xj). Then can be denoted simply as pi without 
any confusion. There are no moduli of holomorphic bundles over the sphere, therefore, we 



a 



a 



where we have used Riemann-Roch theorem in the form 



dim(M a ) - dim(H a ) = N 2 (g(E a ) - 1), 
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can assume that d" A is just the d operator in the line bundle C of the degree deg(E). The 
moment map condition: 

i 

is easily solved: 

4>(x) = -== - — — — - 

% 

Notice, however, that every holomorphic bundle over P 1 has an automorphism group 
GLjv(C), which acts nontrivially on <j> as well as on p^. In fact, the reduction with respect 
to this subgroup is forced by our equation: the sum of residues of cf> must vanish, giving 
rise to the constraint: 

i 

which is nothing but the moment map for the GLn(<E) action. 
Our Hamiltonians in this case boil down to 

Hi, a , b = Res Xa (x - Xa^Trtfix), 

where 1 < 6 < /, 1 < a < L. These Hamiltonians (essentially #2, a, 2) are called Gaudin 
Hamiltonians. 

H haA = Tr(p a ),H 2 , aA =Tr(p 2 a ), 
Tr (p b p a ) 



H 2 ,a,2 = 



etc. 

Spin Calogero-Moser and Rational Ruijsenaars Systems 
Now consider genus zero component E a with double points. Let us decompose the 
set of marked points X a as 

X a = SUTUa(T) 

where t G T and a(t) G cr(T) are mapped to x£a f\ while the restriction of tt onS is 
surjective. We denote pLa = Pu daa^ = 9t- We have: p^jh} 1 = —Ad*(g t )p t . Solving the 
moment map condition as before, we get: 



p t Ad*{g t )p t 



(x-x s ) ^(x-xt) (X-X ff ( t )) 



Sum of the residues vanishes: 



s t 

The group GL N (€) of automorphismes acts on the data p s ,Pt, 9t as follows: 

9t -> 9~ X 9t9,Vt -> 9~ X Vt9,Vs -> £~Vs£ 
Now let us specialize to the case, when #T = 1. The moment map condition will give 

p t - Ad*(g t )pt = ~^2ps 

s 

We have two options: 

we can parametrize the quotient by the action of the group GLn(<C) either by fixing 
the conjugacy class of g t , or the one of p t . 

Consider the first option. Generically one can diagonalize gt, and there will be left a 
group of diagonal matrices, which will act nontrivially on p t and p s 's. Let us denote the 
eigenvalues of gt by 

gt ~diag(e z \...,e ZN ) 
Then in the basis, where g t is diagonal, p t has a form: 

with the further condition 

^2(ps)u = o 

s 

for any i. This condition has an elliptic nature, as we will see, and has a very natural 
origin: the double point on the sphere comes from the pinching the handle. 
Explicitely calculating #2,4,2 we get: 



13 



sink 2 ( i 3 ) 
with 

Sij = 'y^ j (Ps)ij(Ps')ji 



This Hamiltonian describes the particles with spin interaction. The way spins come 
out will be clear later. Right now we simply claim that Sij is the relevant spin interaction. 

Now let us investigate another option - namely, we diagonalize p t . For simplicity we 
shall treat the case #S = 1. We have: 

p t = diag(0i, ...,6 N ) 

and 

(9t)ij(0i-0j) = (Ps)ij 

where p s = gtp s - Now we make a further restriction. Suppose, that for some v e C* 
p s — Idu has rank one. 
Then, 

p s = uld + u <E> v 

where v G (<C N )*,ue<E N . 

Therefore, we can solve for g t : 

[9t)ij ~ Q ._ Q ._ V 
and Ui = (gtu)i. Thus, we obtain a linear equation: 

UjV 

, t - Oj 



E a 3 u 3 = i 

Oi -6t-v 

for any i. The solution is 



P(0i + i/) 

UnVa = 7- — 

where TP (9) = U(°-°i)- 

Finally we introduce the coordinates Zi, defined by 

e Zi Ui = Ui 

The Hamiltonians we can consider in this approach are the characters of g t . We have: 

E «" II 

/C{l,...,iV},#/=fc 1 3 

The system we get is called Rational Macdonald system. If, instead of taking G = 
GLn(€) we would consider SU(N), we will get what is caled a rational Ruijsenaars model 
[R], [RS], [NG2]. 



10 



4-2. Elliptic models 

The next interesting example is the genus one component. Again we omit label a and 
Pap S e ^ s re pl ace d by pi. Generic holomorphic bundles over the torus are decomposable 
into the direct sum of the line bundles: 

E = ®iLi£i 

Therefore, the moduli space M a can be identified with the AT'th power of the Jacobian 
of the curve, divided by the permutation group action. Let us introduce the coordinates 
zi, . . . ,zn on M a - They are defined up to the elliptic affine Weyl group action. Let r be 
the modular parameter of the elliptic curve. Then there are shifts 

Z a —>Z a + 27TV-1 ~ 

T — T 

with m a ,n a G Z, induced by the gauge transformations 

(a- fn r-r n a( x - x ) + m a (xf - xt) 

T — T 

as well as permutations of z^s. Up to these equivalences z^s are the honest coodinates. 

No double points 

First, we dispose of the case, when n doesn't map two points to one. Now the moment 
map condition is 

dfaj + (zi - Zj)4>ij + ^2(Pa)ijS 2 (x a ) = 

a 

The solution of this equation yields a Lax operator 0: for a ^ b 

_ exp(^ J f5f ) a ( Zij + x -x a ) 



13 ~ 27T, 



v^T ^ {PaUj a{ Zij )a{x - x a ) 
where we have denoted for brevity: Zij = — Zj. For the diagonal components we get 

(pii(x) =Wi + ^2(Pa)ii((x ~ Xa) 

a 

In these formulas a and £ are Weierstrass elliptic functions for the curve with periods 1 
and r. The condition of vanishing residues sum is 

a 
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for any a = 1, 



. , N. Now we can compute our Hamiltonians. We have: 



47r 2 Tr(/) 2 (x) = 

Expanding this expression as: 

4ir 2 Tr(j) 2 (x) = p(x - x a )H 2>2>a + ((x - x a )H 2Xa ) + H 2jQ 

a 

we obtain: 

H 2Xa = Ttp 2 a 

as it was expected, 

H 2 ,l, a = ^2w i (p a ) li + ^ (Pa)ii(Pb)iiC(%a ~ %b) + 
i b=£a;i 

These Hamiltonians we shall name as Elliptic Gaudin Hamiltonians. 
The next interesting one is: 

#2,0 = ^2 W i + ^ (Pa)ij(Pa)jip(Zi ~ Zj) + ^ (Pa)ti(Pb)ii(p(x a ~ X b ) ~ ( 2 (x a - X b )) + 
i i=£j',a i;a=£b 

£ ^ (ga " at) (Pa)^(Pt)i ^ r ( Xrr a - X r\ iCiXa ~ Xb + ^ ~ C( ^ )} 
. , . ,, <->\ z ij) (J \^a Xb) 

i^j;a=lb 

Higher Hamiltonians provide us with the integrals of motion of this model. 

Double points 

In the case, when there are the double points, the formulas for the Lax operator and 
Hamiltonians are nearly the same, the only difference is in the conditions on the p a 's. 
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4-3. Spins and coadjoint orbits 

In this section we shall map the notations p a for the Lie algebra elements to the spin 
notations Si, which were introduced in the begining of the paper. 

No double points 

For simplicity, we first consider the case without double points on the component C a . 
First of all, we shall fix the conjugacy classes of all p a 's, that is we restrict ourselves onto 
the invariant subvariety: 

TrPa = Pa,n = COUSt 

Therefore, each p a will represent a point on the coadjoint orbit O a of SL N (€). Again for 
simplicity we assume, that this coadjoint orbit is of the generic type. One can represent it 
as follows. 

Fix a point x a and denote p a simply as p. Introduce a sequence of vector spaces 

S r c . . . c S° 

and consider the space of operators 

Ui : S l ^S l+ \V t : S l+1 
with the canonical symplectic form: 

i 

This form is invariant under the action of the group 

H=*UGL{£ i ) 

by the changes of bases. Therefore, one can make a Hamiltonian reduction at some central 
level of the moment map. Formally it amounts to imposing the constraints: 

U r -iV r -i = Crldg^ 

U i - 1 V i - 1 -V i U i = CiId £i 

for i = 2, . . . , r — 1. Here the complex numbers d are related to the eigenvalues of the 
matrix p via: 

Ai — Xi-i = d 
13 



the multiplicity of Aj equals dim£ l — dim£ t+1 . Finally, p = V U . 

now for each point x a we have a pair of matrices U a ,V a , such that p a = V a U a . 
Then 

(Si) ab = (u a )i ® (v b y e el ® 

Therefore, in our formulas for Hamiltonians we could replace 

(Po)y(Pa)j* -»• Tr £l S i S 3 

Of course, for the products p a pb there is no such interpretation. Thus, Si operators 
in the Gaudin model are just the matrices p a . 

5. Action-Angle Variables 

We first recall the construction of Hitchin in the case of the compact curve of genus 
g>2. 

Given a point in the moduli space of Higgs bundles one can construct a spectral curve 
ScP(T*E©l): 

TZ(x, A) = det((p(x) - A) 

where A is a linear coordinate on the fiber of the cotangent bundle T*E. This curve is 
well-defined, since the equation, which defines it, is gauge invariant. 

The curve S is iV-sheeted ramified covering of E, its genus can be computed by the 
adjunction formula or using Riemann-Gurwitz theorem. 

^)=iV 2 (^-l) + l 

which agrees with the dimension of the moduli space of stable bundles over E. Denote by 
p the projection S — > E. 

Given a stable bundle £ over E we can pull it back onto S. There is a line subbbundle 
C C p*£, whose fiber at generic point (x, A) is an eigenspace of <j>(x) with the eigenvalue 
A. Conversely, given a line bundle £ on E, on can take its direct image, which (again at 
generic point) is defined as 

£x = ®yEp~ 1 (x)£"y 

Therefore, under the flow, generated by the Hitchin Hamiltonians, the C changes and 
it can be shown, that these flows extend to the linear commuting vector fields on the 
Jacobian Jac(S) of S. 
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Thus, the linear coordinates on Jac(S) are the coordinates of the angle- type, whereas 
the intregrals of A over the corresponding cycles in S give the action variables. 

The relevance of this construction in CFT still waits to be uncovered. Presumably, 
it corresponds to Knizhnik's idea [K] of replacing the correlators of the analytic fields on 
the covering of the Riemann surface, by the correlators on the underlying Riemann surface 
with the insertions of additional vertex operators. 

Let us also remark, that quite analagous construction was invented by Krichever in 
in [Kr] in connection to the elliptic Calogero-Moser System. 

5.1. Degeneration of the spectral curve 

We shall adopt the same definition of the spectral curve in the case of degenerate E. 

Obviously, the normalization of S can be also decomposed as the disjoint union of the 
components S a , labeled as the components S a and S a covers E a with some fixed branching 
at the points x l a . Indeed, the behavior of 4> a near the point x % a is known, since the residue 
is known. Let us fix the conjugacy class of p^g. Suppose, that it has k\ eigenvalues of 
multiplicity 1, k l 2 eigenvalkues of multiplicity 2, and so on. Since near the point x l a <p a 
behaves like: 

ij 

(p a {x) ' 



x — x l a 



for appropriate j and /?, then A behaves like 

Pm 



X m (x) 



x - x\ 



where p m is the m'th eigenvalue of p^. Following [BBKT], we can find 

oo 



,1 



m—1 



points Pi, ... , P N ij above x l a , such that the local parameteres Zi, . . . , Z N ij near them are 

aft af3 

defined: 

Z m = A m (x) ^- — - 

The discriminant A a of (p a is a meromorphic iV(iV — 1) -differential on E Q . At each 
point x l a it has a pole of the order 
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The zeroes of A Q dertermine the branching points of the covering 

S a — > E Q 

The number of the branching points equals, therefore, to 

2N(N-l)(g(E a )-l)+J2ol t 

i 

The genus of S a can be computed with the help of Riemann-Hurwitz formula, which gives: 

g(s a ) = i + N*( g -i) + ±J2°« 

i 

Now the Hamiltonian flow due to our Hamiltonians produces a motion of the line 
bundle over S a and it covers the Jacobian of the completed curve S a , therefore, the 
coordinates of the particles will be determined by the same equation: 

i 

as in the simplest one-punctured case. Here 9 is a t/ieta-function on the Jacobian of S a , 
and U defines an imbedding of the moduli space of the holomprphic bundles over E a into 
the Jacobian, as we have described it. 

The details of the recontsruction of all angle-type variables will be published elsewhere 
[NG3]. Remark, that this problem was solved for one-punctured elliptic curve for specific 
orbit in [BBKT]. 

6. Formulas for general case - genus g curve with L punctures 

In this section we consider only one component E of the stable curve. We assume that 
it has genus g and L punctures. We also assume, that E has no double points. 

Using the formulas for the twisted meromorphic forms on the curve, quoted in [I] , we 
can easily write down the formula for the solution of the main equation 

K,0]+^M 2 (^) = O 

i 

In order to do this, we choose the following coordinatization of the moduli space 
J\f of holomorphic bundles over E. Namely, over the open dense subset of Af one can 
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parameterize the holomorphic bundle by choosing a set of g twists: elements of the complex 
group G, assigned to the A - cycles of S. More presicely, let us fix the representatives a k , 
k = 1, . . . , <7, of the A-cycles and let E be the surface E with the small neighborhoods of 
ak removed. Topologically E is a sphere with 2g holes. 

The boundary of the neighborhood of ak consists of two circles a^. In order to 
glue back the surface E one has to attach the projective transformations 7^, which map 
a k~* a k~- These transformations generate Shottky group. On the sphere one can find 
such a gauge transformation h, that 

d" A = hr x dh 

Obviously, 

K9k{x))\- = H k (x)h(x)\ a + 

k k 

where Hk is a holomorphic G -valued function, defined in the vicinity of a~^. Generically 
one can find constant representative of Hk (this is a Riemann-Hilbert problem). 

Once such a gauge h transformation is chosen, the equation for <p can be restated 
in words as the following: find a meromorphic form on E, which satisfy the following 
requirements: 

in the vicinity of x^: d> ~ Pi 

twisting: 4>(^k(x))d^k(x) = AdH k 4>(x)dx 

The answer can be conviniently written in terms of the Poincare seria ([I]): 
Introduce 

u k [x ] G fi^CP 1 ) <g> End(LieG), 
6[x, x ] G fi^CP 1 ) <g> End(LieG), 

where x, Xq G CP 1 , 

u k [x ] z dz = J2MH- 1 )dlog l{z) - lk{Xo) 

9[x,x } z dz = Y Ad{H- 1 )dlog 1 } Z '~ X 

l{z)-x 

where the sum runs over all elements of the Shottky group (it is a free group with g 
generators - loops around 5-cycles. One assigns to the "word" 
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an element of G: 

In this formulas Xq is an auxilliary point on the sphere. 
Finally, the solution has the form: 

g l 

k=l i=l 

The momenta Wk are defined from the condition: 

/ </> = w k 

Ja k 

Easy computation shows, that the symplectic form becomes the sum of the symplectic 
forms on the orbits, attached to the points Xi and of the g copies of the Liouville form on 
the T*G, where the momentum for H k is w^. 

The final remark concerns the vanishing of the residue at the point xq: 

MH^iwk) - w k + = 

k i 

This equation we met in the degenerate form in the section §4.1. It has the meaning of 
the moment map for the action of G on the product of g copies of T*G and the coadjoint 
orbits of pi. 

7. Applications to quantization 

It is straightforward to quantize our models. When the conjugacy classes of p 1 ^ are 
fixed, their quantum counterparts become simply the generators of the group, acting in 
the correspondning representations of G. The condition on the residues of <p a and <pp at 
the double point gets translated to the fact, that the representations, sitting at the points 
x l a and Xp, belonging to one double point, are dual to each other. 

The pinched handle corresponds to the regular representation of the group, and the 
corresponding generators p l ^ a and p^ a are left- and right- invariant vector fields on the 
group. 

Then the Schrodinger equations for the wavefunctions coincide with the critical level 
Knizhnik-Zamolodchikov-Bernard equations [KZ], [Be], [Lo], [EK1],[EFK]. The result of 
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the quantization should follow from the degeneration of the Beilinson-Drinfeld construction 
[Beil]. 

Also, it would be nice to realize the meaning of the generalized KZ equations of [Ch] 
along the lines of our approach. As far as it seems now, these equations are inspired by 
the occasional fact, that the Hamiltonians we have written for the punctured elliptic curve 
are almost symmetric under the exchange :zi <-> x a . 

Finally, note, that using the results of [I] one can easily write down the quadratic 
Hamiltonians for an arbitrary curve (unfortunately, at the moment only in terms of the 
covering of the open dense subset of the actual phase space) , while [FV] allows one to get 
the expression for the wave- functions of the elliptic Gaudin model in terms of the solutions 
of the Bethe Ansatz-like equations. 

When the paper was completed we have been notified about the recent paper by B. 
Enriquez and V. Rubtsov [ER] on the related subject. We would like to thank authors of 
[ER] for their comments. 



19 



REFERENCES 



[BBKT] O. Babelon, E. Billey, I. Krichever, M. Talon, Spin generalization of the Calogero- 
Moser system and the Matrix KP equation, hepth/9411160 
[Beil] A. Beilinson, V. Drinfeld, Quantization of Hitchin's fibration and Langlands program, 
A. Beilinson's lectures at IAS, fall 1994 
[Be] D. Bernard, Nucl. Phys. B303 (1988)77, B309 (1988), 14 
[C] F. Calogero, J.Math.Phys. 12 (1971) 419; 
[Ch] I.Cherednik, Difference- elliptic operators and root systems, \hep-th/9410188{ 
[EF] P.Etingof, I.Frenkel, Central extensions of current groups in two dimensions, 
hepth/9303047 

[EK1] P.Etingof, Quantum integrable systems and representations of Lie algebras, 
hepth/9311132 

[EK2] P.Etingof, A.Kirillov, Jr., On the affine analogue of Jack's and Macdonald's polyno- 
mials, Yale preprint, 1994 

[EFK] P. Etingof, I.Frenkel, A. Kirillov-Jr., Spherical functions on affine Lie groups, Yale 
preprint, 1994 

[FV] G. Felder, A. Varchenko, Integral formula for the solutions of the elliptic KZB equa- 
tions, hepth/9503001 

[NG1] A. Gorsky, N. Nekrasov, Elliptic Calogero-Moser System from Two Dimensional Cur- 
rent Algebra, hepth/9401021 

[NG2] A. Gorsky, N. Nekrasov, Relativistic Calogero-Moser model as gauged WZW theory, 
Nucl.Phys. B 436 (1995) 582-608 

[NG3] N. Nekrasov, A. Gorsky , Integrable systems on moduli spaces, in preparation 
[FF] B. Feigin , E. Frenkel, N. Reshetikhin, Gaudin Model, Critical Level and Bethe Ansatz, 
CMP 166 (1995), 21-62 

[G] M. Gaudin, Jour. Physique, 37 (1976), 1087-1098 

[H] N. Hitchin, Duke Math. Journal, Vol. 54, No. 1 (1987) 

[I] D. Ivanov, Knizhnik-Zamolodchikov- Bernard equations on Riemann surfaces., |hep 



th/9410091 



[KKS] D.Kazhdan, B.Kostant and S.Sternberg, Comm. on Pure and Appl. Math., vol. 
XXXI, (1978), 481-507 
[KZ] V. Knizhnik, A. Zamolodchikov, Nucl. Phys. B247 (1984) 83 
[K] V. Knizhnik, Analytic fields on Riemann surfaces, CMP 112, (1987), 567 



20 



[Lo] A. Losev, Coset construction and Bernard Equations, CERN-TH. 6215/91 

[OP] M.Olshanetsky, A.Perelomov, Phys. Rep. 71 (1981), 313 

[M] J. Moser, Adv.Math. 16 (1975), 197-220; 

[Kr] I. Krichever, Funk. Anal, and Appl., 12 (1978), 1, 76-78; 14 (1980), 282-290 

[RS] S.N.M. Ruijsenaars, H. Schneider, Ann. of Physics 170 (1986), 370-405 

[R] S. Ruijsenaars, CMP 110 (1987), 191-213 

[S] B. Sutherland, Phys. Rev. A5 (1972), 1372-1376; 

[V] E. Verlinde, Nucl. Phys. B 300 (1988) 360 

[ER] B. Enriquez, V. Rubtsov, Hitchin systems, higher Gaudin operators and r-matrices, 
|alK-geom/950301C 



21 



